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We present a computation of the decay constant Fb^ in quenched QCD. Our strategy is to combine new precise 
data from the static approximation with an interpolation of the decay constant around the charm quark mass 
region. This computation is the first step in demonstrating the feasability of a strategy for Fb in full QCD. The 
continuum limits in the static theory and at finite mass are taken separately and will be further improved. 



1. Introduction 

A non-perturbative precise computation of Fb 
in full lattice QCD with reliable errorbars would 
be a major achievement with important phe- 
nomenological implications for CKM physics and 
the search for new fundamental processes. Such 
a computation has to meet three main obsta- 
cles: chiral extrapolation, the heavy b-quark and 
the large computational cost due to unquenching. 
We avoid chiral extrapolations for the time being 
by setting the light quark mass to the strange 
mass P , thus addressing Fb^ , which is of interest 
in itself, for example in Bg — Bg mixing. A strat- 
egy to deal with the second problem has been sug- 
gested in |2I3| . The idea is to use HQET and the 
1/m expansion after a non-perturbative match- 
ing to QCD in finite volume. No large lattices 
are required and the method can thus be applied 
to quenched as well as full QCD. This strategy 
seems to be viable now thanks to the consider- 
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able improvement of the statistical precision for 
computations in the static approximation to QCD 
that has been found in (4t5) The required order 
of the l/m expansion is unknown. However, we 
hope that the size estimates of the terms in the 
HQET presented in W. are correct and the lin- 
ear term proves sufficient to keep the systematics 
under control. This can with the present com- 
putational power only be tested in the quenched 
approximation. 

Here we compare the renormalization group 
invariant matrix element $f^Qj of the static ax- 
ial current to relativistic data around the charm 
quark mass. Later we will also compute the l/m 
term in the static approximation to complete the 
validation of our strategy for Fb^ ■ We use an in- 
terpolation of our relativistic data and of $rqi 
to obtain a precise quenched value of Fb^ in the 
continuum limit. Previous results for Fb^ have 
been reviewd in [7] and a new computation has 
recently appeared [S]. 



2. Numerical Results 



The renormalization group invariant matrix el- 
ement ^^Qi at infinite mass is related to the pseu- 



doscalar decay constant Fpg at finite mass by a 
matching factor Cpg, 
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Cps can be expressed as a function of the renor- 
malization group invariant heavy quark mass M , 

M 

X l/log(-r ) < 0.62 => 



x--{l - 0.06814a;- 0.08652 

+0.07939 x^). (2) 

This expression is a simple and accurate 
parametrization of Cps obtained by integration 
of the entering renormahzation group functions as 
explained in W, where the anomalous dimension 
of the static axial current is taken to 3 loops ^U] ■ 
Its uncertainty is estimated to be smaller than 
two percent which is half of the difference between 
the 2-loop and the 3- loop result. 

Our computational setup has been explained 
in |11I12| . In particular we use Schrodinger 
functional boundary conditions and employ non- 
perturbative 0(a) improvement. 

We obtain the decay constant at five meson 
masses from « 1.7 GeV to « 2.6 GeV at four lat- 
tice spacings a ranging from 0.1 fm to 0.05 fm. 
Since the data at the different lattice spacings 
are not at exactly the same pseudoscalar masses, 
we interpolate them linearly in l/rnpg and eval- 
uate the interpolating functions at a few meson 
masses. These numbers can then be extrapolated 
linearly in to the continuum limit. This has 
to be done with care however, since the lattice 
artifacts depend on the mass. To avoid large 
discretization errors in the slope of the contin- 
uum extrapolation we thus follow the experience 
from perturbation theory 13 and leave out more 
and more coarse lattices at larger masses. At the 
largest masses only the finest lattice can be used 
and we take that point as the continuum limit. 
We add to its error the difference to the contin- 
uum value obtained from a linear fit from the two 
finest lattices. The continuum limits for two pseu- 
doscalar masses are illustrated in figure which 
shows the case with the two coarsest lattices left 
out in the upper part and the case where only one 
lattice contributes in the lower part. 
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Figure 1. Continuum extrapolation of the decay 
constant Fps^/mpg at two meson masses from 
two resp. one lattice. 

Our result is still preliminary since we will add 
an even finer lattice with a w 0.03 fm (cmp. \V2\) 
and since the details of the continuum extrapola- 
tions are still being discussed. Finally we choose 
five points close to our original data to be used 
in the interpolation between the static and the 
relativistic case. 

To compare with we still have to compute 
the renormahzation group invariant masses M 
that are needed for the evaluation of the match- 
ing factor Cps. Here we follow exactly ^3] for 
meson masses in the range considered. 

We take the static result from 



<^t^\ = 1.74(13), 



(3) 



which has been obtained in the continuum 
limit from three lattice resolutions with a ~ 
0.1 fm ... 0.07 fm with a new static action that 
uses HYP smeared links |15| . Furthermore, wave 
functions for the states at the Schrodinger func- 
tional boundary have been used together with 
an elaborate technique to reduce the contribu- 
tion from the first excited state ^ and the 
renormahzation factor has been computed non- 
perturbatively ^ . 

Our results are shown in figure |5] as a func- 
tion of the inverse pseudoscalar mass. At those 
points where the number of lattices taken into 
account in the continuum extrapolation changes, 
there is a small systematic effect which explains 
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Figure 2. Comparison of the renormalization 
group invariant matrix element of the static axial 
current with relativistic data around the charm 
quark mass. 

the zigzag behaviour of the 1-sigma band around 
the relativistic data. If we combine the five points 
selected above and the static point by an inter- 
polation in 1 /mps we can obtain F^^ . To this 
end we use Mb = 6.96(18) GeV P] and A^=° = 
238(19) MeV to evaluate the matching factor 
Cps at the b-scale. Both, a linear and a quadratic 
interpolation lead to 

Fb, = 206(10) MeV. (4) 

Here we have used me^ — 5.4 GeV and set the 
scale with vq = 0.5 fm. 

Furthermore we notice, that the data are well 
described by the linear interpolation that is dis- 
played in the plot. The slope of the fit is —2.1(5) 
while without the constraint through the static 
approximation we would obtain —2.1(1.1). To be 
able to compare this result with the slope com- 
puted from the HQET it is desirable to get this 
slope even more precise. 

As in jJT] we can use the slope of the interpo- 
lation to get an estimate of the quenched scale 
ambiguity of Fb^, i.e. of the effect of changing rg 
within 10%. Under this change Fb^ changes by 
12%. Note that the true quenching error can of 
course only be estimated by an unquenched cal- 
culation. 

3. Conclusion 

Our comparison of the static and the relativis- 
tic results shows that at the current level of statis- 
tical precision nothing contradicts the hope that 
the 1/m expansion including only the linear term 



is enough to compute Fb^ from the HQET in full 
QCD at a sufficient precision. To further validate 
this strategy we will compute the 1/m term in 
the quenched approximation from HQET. 

The combination of with our data at 

finite mass results in a new value for Fb^ in 
quenched QCD, Fb, = 206(10) MeV, including 
all the systematic errors up to quenching. The 
quenched scale ambiguity is estimated by 12%. 
Since we interpolate our data, this is a stand 
alone result that is almost independent of any 
effective theory, as is indicated by the fact that 
a linear and a quadratic function yield identical 
values. 
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